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Abstract 

We have introduced a spin-isospin dependent three-dimensional approach for formulation of 



C^ I the three-nucleon scattering. Faddeev equation is expressed in terms of vector Jacobi momenta 

and spin-isospin quantum numbers of each nucleon. Our formalism is based on connecting the 
transition amplitude T to momentum-helicity representations of the two-body t-matrix and the 
deuteron wave function. Finally the expressions for nucleon-deuteron elastic scattering and full 
breakup process amplitudes are presented. 
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I. INTRODUCTION 

During the past years, the three-dimensional (3D) approach has been developed for few- 
body bound and scattering problems [1-15]. The motivation for developing this approach 
is introducing a direct solution of the integral equations avoiding the very involved angular 
momentum algebra occurring for the permutations, transformations and especially for the 
three-body forces. Furthermore above the pion production threshold and the GeV region the 
three-nucleon (3N) scattering with the 3D formalism is simpler for more complex few body 
system by providing a strictly finite number of coupled 3D integral equation. Conceptually 
the 3D formalism considers all partial wave channels automatically. 

In the case of the 3N scattering the Faddeev equation has been formulated for three 
identical bosons as a function of vector Jacobi momenta, with the specific stress upon the 
magnitudes of the momenta and the angle between them pj, [2|. Adding the spin-isospin to 
the 3D formalism is a major additional task, which will increase more degrees of freedom. 
This carried out by Fachruddin et al. by considering leading order of Faddeev equation for 
breakup process. Recently a 3D formalism based on operator form has been introduced for 



formulation of the 3N scattering [16|. In this paper we have attempted to formulate the 
full breakup and the elastic scattering of nucleon-deuteron (Nd) by including directly the 
spin-isospin degrees of freedom. To this end we have formulated the Faddeev equation for 
the 3N scattering with the advantage of using the helicity representation of NN forces such 
as the AV18, Bonn-B and the Chiral potentials j^- jlTj. 

This manuscript is organized as follows. In sect. Hi] we have derived Faddeev equation in 
a realistic 3D scheme as a function of Jacobi momenta vectors and the spin-isospin quantum 
numbers. In sect. IIIII we have derived two expressions for amplitudes of elastic scattering 
and breakup process respectively. Finally in sect. llVt a summary and an outlook have been 
presented. 

II. FADDEEV EQUATION FOR THE 3N SCATTERING IN THE 3D APPROACH 

Faddeev equation for the three identical particle is given by [l8|: 

T\ct>) = tP\4>) + tGoPT\4>), (1) 



where P = P12-P23 + P13P23 is the sum of cychc and anti-cychc permutations of the three 
nucleons, t denotes the NN transition matrix determined by a two-body Lippman-Schwinger 
equation and Gq is the free 3N propagator which is defined as: 
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(2) 



where Hq and E are the free 3N hamiltonian and the total energy in the center of mass 
frame respectively, Ed is the deuteron bounding energy and qo is the relative momentum of 
the projectile nucleon to the deuteron. In order to solve eq. O in the momentum space we 



introduce the 3N free basis states in a 3D formalism as 



d 



m: 



|pq7) = |pq m^^m^2 "^,3 mij mi^mig ) = Iqm^^mi^) Ipm.^m.gmi^mig), (3) 

The basis states involve two standard Jacobi momenta p and q which are the relative mo- 
mentum in the subsystem and the momentum of the spectator with respect to the subsystem 
respectively [19']. I7) = jm^-^ 771^2 msgrni^m^jmig) is the spin-isospin parts of the basis states 
where the quantities m^- {mt^) are the projections of the spin (isospin) of each three nucleons 
along its quantization axis. The introduced basis states are completed and normalized as: 

Y, dp dq, |pq7) (pq7| = 1, (p'q'Y IpQT) = '^(p' - P)5(q' - q) ^7- (4) 

Projecting the Faddeev equation on to Jacobi momenta leads to: 

(pq7|r|qomOm?vl>f<*) 

= (pq7 1 tP + tGoPT I qomX^f'') 

= (pq7|tP|qom°m°vI/f ) + (pq7|tG'oPT|qomOm°vl/f^), (5) 

where: 

10) ^ IqomX^^) = |qom°m°)|v^f''), (6) 

is the initial state where m^ (m^) is spin (isospin) projection of projectile nucleon along quan- 
tization axis and \^ d"^) ^^ ^^^ antisymmetrized deuteron state with Md bing the projection 
of total angular momentum along quantization axis. We start by inserting the completeness 



relation twice into the first term of eq. ([5]) : 

(pq7|tP|qom°m°vl/f) 

= E / ^P / rfq'(pq7|t|p'q'7') E / dp" I rfq"(p'q'7'l^|p"q"7") 

7' 7" 

x(p"q"7"|qom°m?vI;f). (7) 

Tlie matrix elements of the permutation operator P are evaluated as [Tj : 

(p"q"7"|^|p'qY) 

+ (5(p" + -q" + q') 6{p' - q" - -C^)6m'' m'Sm''m'Sm''m' ^m'l m' ^m'f m' ^m'f m' ■ (8) 

We have used these relations for the matrix elements of the two-body t-matrix and the 
deuteron wave function: 

(pq7|t|p'qY) = (P"^S2'^S3"^t2"^t3l^(^)|P'K2"^'.3K2"^t3)'^(q~q')'^m,,m^^5^,,mJ^, (9) 

(p"q"7"|qom>°vl/f^) = (p"m;>';m;>;;|vl/f ) 5(q" - qo) 5^,, „o <5„,, ^o, (10) 

where the two-body subsystem energy in a 3N system is ^ = E — -^q^. Substituting eqs. (|H]), 
da]) and (UnD into eq. (CD yields: 

(pq7| tP|qomgm°\E'^ '*) 

= X^ (P?«s2"^s3'^t2'^t3l^(^)l -7rm',2m°mj2m°)(7r'msim',2mt,m^J^f'*) 
K2K2 
+ X^ (P"^s2?^s3"^t2"^t3l^(^)k"^Vs3"^tK3)(-7rW,3m,,mJ3miJ^f^) 

<3™'3 

= ^ < (pms2ms3mi2mi3|t(2;)P23|7rm°m^m°mJ)(7r'msim^mi^mj|\E'^'^'*) 

-|-(pms2ms3mt2rn(3|t(z)|7rm°m'^rn°mJ)(7r'ms-^rn'^mijmJ|P23"'^|^d '') 
= - 5Z (P"^s2"^^3"^t2"^t3l^(2^)(l - P2z)\T^mlm'^m^^m[){T^'ms^m'^mt^m[\-^^''). (11) 

In the last equality we have used the antisymmetry of the deuteron state 1^1'^'') and also we 
have considered: 

1/1 

7i" = 2q + qo> 7r=q+2qo' (12) 



The antisyininetrized two-body t-matrix is given by ^: 

(13) 
where |pms^ms2W-timj2)a is the antisymmetrized two-body state which is defined as: 

\pms^ms^mt-,mt^)a = —^{l - Pi2)\pms^ms:,mt^mt^) . (14) 

Therefore eq. (11) can be written as: 
(pq7|tF|qom°m°^f'*) 
= -^ a(pm^2ms3mt2mj3|t(2;)|7rm°m'^m°mj)a (Tr'm^jm^mj^m^lvl/^-^''). (15) 

Now by inserting the completeness relation twice into the second term of eq. (|5]) and using 
eqs. ([8]) and ([9]) we have obtained: 

(pq7|tGoPr|qom°mOv|>f<^) 

= E / ^P" / V(pq7|t(^)G'o|p"q"7") E / ^P' / ^q (p"q"7"l^|P q'7') 
7" 7' 

x(p'qY|T|qomV?^f^) 

/ / / I J -C/ "T — h ^£ 

m' m' rrij m,j m 4m 

, , I , J E - - — '^ +ie 

m'm'm[ m[ m Am 

X (vi-'q'm'^^ rris-^ m'^g m^^ ruij m^ J T I qom°m°\E'^^'' ) 
, ,, J E-^-^ + ie\ 

m' m'^m^ m[ m 4m V> 

x(7r'qW^^m5^m!,mJ^mi^mJ|P23^r|qom°m°^f'') + {pms^ms^mt^mt^\t{z)P2-i\Tirn'^m'^m[m[) 
= - Y] / d<^- p §3 X {pms2ms,mt,mt,\t{z){l - P23)\nm m'^m[^m[) 

m'm'gmim]. m 4m 

v,/~' ' ' ' ' '\'-n\ 0,TcAfrf\ 

E / t/q'- p ^^2 -a{pms2ms,mt2mtMz)\^K^K^'h^t)a 

n'^m^ m[ m Am 



where we have considered: 

7^=2^ + ^' 7r' = q+-q'. (17) 

Final expression for the Faddeev equation is exphcitly written as: 

= - ^ <^ a{pms2ms^mt2mt^\t{z)\nm'lm'^m^m[)a{n'ms^m'^mt^m[\'^^'') 



(18) 



(7r'q'm'^^msimgm^^mj^m^|T|qom°m°\l/^ '*) 

^ ^-;^(g' + q-q' + g'') + ^£ 



Since the transition operator T is needed for all of q values and We know that the two-body 
interaction supports a bound state, which is characterized by a pole in the two-body t-matrix 
at the two-body binding energy Ed-, thus we need to consider this pole, which is located at 
z = Ed- The residue at the pole can be explicitly extracted by defining: 

t-^ = {z- Ed)t\ (19) 

Since the pole of t" will be present in T, we define: 

f={z- Ed)T, (20) 

Therefore eq. (18) can be written as: 

(pqm^jms2m,3mtimt2mi3|f|qom°m°^f'') 
= - 5Z 1 a(p"is2m^3mt2mt3|t(z)|7rmV5"T'Vt)a(7r'm^im'^mtim^|^f'') 

, a{vrns2ms3mt^mt^\i{z)\T(m'^^m'^m[^rn'^)a 



+ E /^q' 



™si™ii 



^-^(?' + q-q' + ?") + 



le 



(7f'q'm'^^ms-^m'^mJ^m(^mJ|T|qom°m°\l'^ '') 
E-^q'^-Ed + ie 



x """^""^^"^-'7-^r : : " ^ "^ ' )■ (21) 



As a simplification we rewrite eq. (21) as: 






\ ^ J Z m'im''rr£m't /__ _. .A vtiM, 






m^rrq ,Mj^ 



V"^ I 1 It f-a ms2ms3mt2mt3 l^P, TT, 2j ra'lms^m'^ra'lmt-^m[\ '^ i MU; I /r)r)\ 

For solving this integral equation one needs the matrix elements of the deuteron wave func- 
tion and the antisymmetrized two-body t-matrix. We have connected these to their helicity 
representations in appendices A and B respectively. 

III. DERIVATION OF THE 3D EXPRESSIONS FOR THE 3N ELASTIC AND 
BREAKUP SCATTERING AMPLITUDES 

A. Elastic scattering 

In the Faddeev scheme the operator U for the Nd elastic scattering is defined by: 

U = PGq^ + PT, (23) 

Differential cross section for Nd elastic scattering in the center of mass frame is given by: 
^ = {2n)\^f^- y |[/^-^^. o(q,qo)P, (24) 

where U^"^^ mOmo(^'^o) ^^ ^^^ elastic scattering amplitude which is defined as: 

U''-''\ o(q,qo) 
= (qm^mi^f''|[/|qom°m°^^'') 

= (qm.mi^f^lPG^^IqomV^^f') + (qm,mi^f^|PT|qom°m°^f''), (25) 

where: 

|qom°m°^f^) = |qomV?)l^f^), qm,mi^f^| = (qm,mi|(^f^|, (26) 



are the initial and the final states respectively. For derivation of an expression for the elastic 
scattering amplitude we start by inserting the completeness relation twice into the first term 
of eq. 1^5^ as: 

p'2 3q'\ 



x(p"q'YI^|p'qY)(p'qY|qom°m?vi/f^). (27) 

with considering eq. ([H]) and following relations: 



(p'q'7lqom°m?vl;f^) = {p'm',y,y,y,J<I^f^)5{q! - qo) S^,^^o6^,^^o, (29) 
Equation ([27]) can be written as: 

(qm^mt^^'' I PGt^'^ I qom°m°^^'' ) 

K2K2 

<3"'t3 

= -2{E (g2 + q ■ qo + go)) ^ {^''f''\Timlm^mlm[)(ji'msm',mtm[\^^'''). (30) 



m 



Now by inserting the completeness relation twice into the second term of eq. ( !25|) and using 
eqs. ([8]) and fl28|) we have obtained: 

= E / dp" I rfq"(q"^."^*^f 1p"q"7") T^J^p' J rfq'(p"q"7"l^|p'qV) 



'y" 

X (p q 7 |r|qom^m"W^ ") 



\ "* I J I) /iTc*^dl~ I I I l\ l~l II II l\ n — IrT-il 0,TrA^rf\ 

+ (\l/^ |i23|vrmg^m^m(^mj)(7r qm^^msm^mj^mimji Iqom^m^W^ ) 

^ -2 j: /..'(.f |,K.™>;.™'.>(*'.'™>.™>;.-n,KiTi<,o»»<*-> 

= —2 > / dq (W , " 7rm„ m„m. m.) i .-^^ •(31) 

Final expression for the matrix elements of elastic scattering amplitude has been obtained 
as: 

(qm^mt^f '^ |t/|qom°m°^^^'*) 



+ > dqiW.'^l'Km^^m^'m.m.) ^-^^ P- 32 

A^ J ^\ d \ s, s t, t/ E-^q'^-Ed + ie j ^ ^ 



Tcmlm'gm^'m^) {n'msin'^mtm^ \ \1'^ '^ ) 

0^0,Tf^d\ 



As a simplification we rewrite eq. f p2|) as: 

u " o(q,qo) 






-K) 



\mflrrfi,M. 



x^>''^ , ,(7r')+ y [dc^'^r^ „ I in) <--<Kw^^^'q'^°u . (33) 

m'Jm'J ^'^ y 






B. Full breakup process 

In this stage we have derived an expression for the matrix elements of the breakup process. 
In the Faddeev scheme the operator Uq for the Nd breakup process is given by: 

f/o = (1 + P)T. (34) 

Differential cross section for the Nd breakup process is defined by j:7]: 
da -^ 



MamO-y 

where the matrix elements of the full breakup amplitude Uq '^"^^'^(p, q; qo) is defined as: 

= (pq7|(l + P)r|qomVOvl/f^) + (pq7|r|qomV?^^) 

+ (pq7|(Pi2P23)T|qomX^f') + (pq7l(A3P23)T|qom°m°M/f''), 

(36) 

where |qo?ri°wi^\I'^ '') and (pq7| are initial and final states respectively. By applying the 
permutation operator P12-P23 and P13P23 fo the final state, eq. (l36ll can be written as [Tj: 



^0 IP' qi qoj 



0^0,T,Md\ 



+ ((-2P - 4q)(P ~ 2^''"^"2"^"3"^si'^t2"^t3"^til^|qo"^°"^?^f^) 

13 1 

+ ((-2P + 4q)(-P - ^(^rns^ms^ms^mt^mt^mt^\T\qQmlm^^^!'^'') 

((-|P- iq)(p- |q)m,,m,3ms^mt,mf3mtjT|qomgm^^f'^) 
((-|P + |q)(-P - |q)m,3m,^m,,mt3mf,mt2|t|qomgm°^f'') 



(37) 



10 



As a simplification eq. ( 137|) has been rewritten as: 









+ 



J-m^^m^, m^„mt„mf, mt,A nP+ ^Qi P n* 



^-il^(-p-ip)^-^^ 



4m ' 

(38) 



IV. SUMMARY AND OUTLOOK 

We extended the recently developed formalism for a new treatment of the 3N bound 
state in three dimensions for the Nd scattering [13]. We propose a new representation 
of 3D Faddeev equation for the 3N scattering including the spin and isospin degrees of 
freedom in the momentum space. This formalism stays closely to the bosonic structure 
where the spin and isospin degrees of freedom are ignored. This is an important step forward 
since our 3D formalism avoids the very involved angular momentum algebra occurring for 
the permutations and transformations and it is more efficient and less cumbersome for 
considering the 3N forces. This formalism enables us to realistically handle more complexity 
in 3N scattering calculations. This work provides the necessary formalism for the calculation 



of 3N scattering observables which is under preparation 20 ] 
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Appendix A: Connection of the deuteron wave function to its helicity representation 

In In our formulation, we need the matrix elements of the deuteron wave function 
^flf'' rr> rr, ™ (p) • Wc havc couuectcd these matrix elements to those in the momentum- 
helicity basis. The momentum-helicity basis state which is parity eigenstate and antisym- 



11 



fl: 



metrized is given by j^: 

|p;p5A;t)™ = -L^l-P,,)\p-pS\)^\t) 

= ^{i-vA-)'^')\p;pSX)^\t), (Al) 

Here S is the total spin and A is the spin projection along relative momentum of two nucleons. 
|t) = \tmt) is the total isospin state of the two nucleons, where t is the total isospin and rrit 
is the isospin projection along its quantization axis, which tells also the total electric charge 
of system. For simplicity rrit is suppressed since electric charge is conserved. Pu is the 
permutation operator which exchanges the two nucleons labels in all spaces i.e. momentum, 
spin and isospin spaces. |p; pS'A)^ is parity eigenstate which is given by: 

\p;pSX)^ = ^{l + rj^P^)\p;pSX), (A2) 

where P^ is parity operator, r].,, = ±1 are the parity eigenvalues and |p;p5'A) is the 
momentum-helicity basis state. The normalization of the momentum-helicity basis state 
is worked out as: 

-'"(p';p'>S'A';t'|p;p5A;t)™ 
= (1 - vA-f^') St'tS,^,vJs's{S{p' - p)6x'x + VnH'Sip' + p)6y^,}, (A3) 

and the completeness relation of this state is defined by: 

J2 [dp IP; PSX; tr \ ^"{p; PSX- t\ = l. (A4) 

Inserting the completeness relation in the momentum helicity basis yields: 



1 r 

= I E / dp\pm,,ms,mt,mt,\p';p'lX-QY'' ^''(p'; p'lA; 0|vi>f^). (A5) 

where we have used the deuteron properties, i.e. 5* = 1, t = and even parity, ^^''^{p') = 
^°(p') p'^X; 0|\E'^ '') is the deuteron wave function component in the momentum-helicity basis. 
The overlap of the momentum-helicity basis state with the state IpmsjmsjTTii^mtj) is given 
by [41: 



{pms-,ms.,mt^mt2\p'; pSX; t)' 

X e-'^^^dl^{e) [5{p - p') + r^J{p + p')], (A6) 



1 11 11 

-(1 - ?7^(-)^+*) C(--; mt^mt^) C(--; m.^m,^ Aq) 
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where C is the Clebsh-Gordan coefficient and rff )^{9) are rotation matrices 2l|. Substituting 
eq. (IA6I) into eq. flASP yields: 



A=-l I 

+e-'^o(-+-)<,(vr-^)$f(-p)|. (A7) 

By considering these relations: 

<a(^ - ^) = (-)"^'''<-aW, <^f (p) = -<^-'^(-p), (A8) 

Equation (]A7|) can be written as: 

A=-l 

(A9) 

The azimuthal dependency of ^^ '' (p) is: 

<l>f <* (p) = e^*^^^<l>f '^ (p, 6') (AlO) 

Finally the connection between the \E'^^ ^^ ^^ ^^ (p) and those in the momentum-helicity 
basis state is given by: 

A=-l 

(All) 



It should be mentioned that the deuteron wave function component ^^"^(pjO) obeys a set 
of coupled equations which are solved numerically in ref . [2^ . 

Appendix B: Appendix B. Anti-symmetrized NN t-matrix and connection to its 
helicity representation 

The connection of Anti-symmetrized two-body t-matrix to those in the momentum- 
helicity basis is given in ref. 7j , here we prepare this connection according to the notation to 
be used in our work and then we have derived an expansion in momentum-helicity basis for 

13 



two-body t-matrix for z — )■ E^. Based on moinentum-helicity basis states the NN t-matrix 
element is defined as: 

t5;f/(p, p';z) ^ -'^(p;p5A;t|t(^)|p';p'5A';t)-^ (Bl) 

As shown in ref. [TJ, the selection of p' parallel to the z-axis allows, together with the 
properties of the potential, that the angular dependencies of the NN t-matrix elements can 
be simplified as: 

= e-^^^vv' e'^''^pp' tf J {p,p', cos epp>;z) 

= e*^"''^-'-^^-') tfj ip,p', cos eppr,z), (B2) 

The direction hppi can be determined by the spherical and the polar angles 9pp' and ippp', 
where 

cos 9pp' = cos 9p cos 9p' + sin 9p sin 9p' cos((/)p — (t)p'), 
sin 9 pp'e^^'^pp' = — cos 9p sin 9p' + sin 9p cos 9pi cos{(f)p — (ppi) + i sin 9pSm{(l)p — cpp'), (B3) 

and the exponential factor e**-^ <t>pp'-^^) is calculated as: 

EL-s Dfjxi<Pp dp 0)D%Wp' Op' 0) 



e'^^pp' 



e^i^'^pp'-^^p') 



Dyxi'Ppp'dpp'O) 

^%_^e^^('^p-'^.')dUOp)df,AOp' 



4a(V) 

(B4) 

In the above expressions, Dff^{(f>p9pQi) are the Wigner D- functions. Finally the connection 

between the t-matrix elements a(p?7T,s;^r?T,s2 ?TT,j^mt2|t(^)|p'm'^j?TT,g2 '^ti'^ta)" ^"^^ those in the 
momentum-helicity basis, namely t\f}{'p,'p'\ z), is given as: 



a{^msmsmt^mt^\i{z)\^'m'^^m'^^m[m[^)a 



xe 

SlTt 



11 11 



E<A(^p)4A'(^p')ilf' (P, P'; z). (B5) 



X 
AA 
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It should be mentioned that the fully off-shell two-body t- matrix t\f}{j),p' , cos 9pp'] z) obeys 
a set of coupled Lippman-Schwinger equations which for S* = it is a single equation and 
for S" = 1 it is a set of two coupled equations which are solved numerically in ref. [? ]. 

In this stage we have derived an expression in the momentum-helicity basis for NN t- 
matrix for z — > E^. It is clear that: 

lim {z - Ed)t = y|^f^)(^f^|V. (B6) 

2 >Ea 

Projecting this equation on momentum-helicity basis states yields: 



110/ 
\X'\ 



]im {z-E,)t\'^,{p,p' 

z — >Ea 



In [dp"'^{p;plX;0\V\p";p"lX";Oy'^ '^{p";p"lX";0\^f^) 
X" -^ 

-J2 [dp"{^f^\p";p"lX";Oy^ ^"(p"; p"lA";0|y|p';p'lA';0)i" 

A" "^ 

\ Y: J dp"v^'Ap, p") <^p{p") 1 1 i E / dp'yy'W^ p') ^i^(p") I • m 

By considering vector q' along z axis we have obtained: 



"■4 

A" 



e'''^ \im(z-E,)tl'^,ip,p',x) 

z — >Ea 



= (lE / w^-v,^^°(p,p")$;:^-(/,x") 

I A" -^ 

X { i E / rfp"e-^^^^"e^^'-"n^4 (/, p', x") <!>^np", x") I , (B8) 

where we have used the azimuthal behavior of the potential and the deuteron wave function 
as: 

^^,;'{p") = e'^'^^"^'^,f{p\x") (B9) 
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Equation (IBSP can be rewritten as: 



z — >Ed 



X" ^ 

X" ^ 
\ E l^ '^P"P'" j\ d^"^'xy'''iP^ ^y^ ^") ^x''iP"^ ^") 

{\^[ ^^"^"' £ '^^" ^^'' °' (^"' ^'' ^") ^'^'' (^"' ^") /" '^^"e^^''"*'''^^"'""^ \ > 



X 

1 4 

A' 



(BIO) 



where we have introduced: 

r.27r 



vfXHp.^.P".^")= / rf^"e-^*^^('^--'VA^?(p,p"), (Bll) 



Finally we have obtained: 



Wm {z-Ed)t]^;^,{p,p\x) 

z — >Es, 



{ i E /°° <^P"P"' /' dx"vf^,;'Hp, x,/, x") ^^,f{p\ x") 

\ "^^'^^^ E [ ^P"P"' /' ^^" ^A''v ip",P', ^") Kf{p", ^") I , (B12) 
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